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Abstract 

In this paper the heat signahng in microtubules (MX) is investi- 
gated. It is argued that for the description of the heat signahng phe- 
nomena in MT, the hyperbohc heat transport (HHT) equation must 
be used. It is shown that HHT is the Klein-Gordon (K-G) equation. 
The general solution for the K-G equation for MT is obtained. For 
the undistorted signal propagation in MT the Heisenberg uncertainty 
principle is formulated and discussed. 
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1 Introduction 



The aim of this paper is to provide a semi- quantum theory of intracellular 
heat transport of cell organelles and vesicles here termed "particles". Nu- 
merous experimental studies have established that this particles are equipped 
with bound motor proteins which move them along microtubules and actin 
filaments For example anterograde transport of particles along micro- 
tubules in nerve axons is mediated by the motor protein kinesin [2]. In this 
system the motion of particles is not continuous but saltatory particles 
are transported for distances of typically ~ 10 nm at a more or less steady 
velocity of ~ l/ims~^ but there are pauses lasting for upward of 1 s in which 
given particle is apparently undergoing Brownian motion and has presumably 
detached from the microtubule or is stuck. 

Recently D. A. Smith and R. M. Simmons |3] developed the reaction - 
diffusion model for the motor assisted transport of intracellular particles. In 
this paper we propose the model for signaling phenomena inside the micro- 
tubules. Considering our results collected in the monograph [5^ we argue 
that the transport phenomena inside the microtubules must be discussed in 
the frame of the quantum transport equations for the diameter of the micro- 
tubules is of the order of nanometers, i.e. the order of the molecules. In the 
paper we develop the Klein-Gordon (K-G) type transport equation and will 
obtain the solution of K-G for Cauchy boundary conditions. 

Recently [TO] the study of the motion of the kinesin motor along mi- 
crotubules using interference total internal reflection microscopy was under- 
taken. This technique bears some similarity to earlier ones designed to study 
the proximal hydrodynamics behavior of polymers to a surface. The laser 
source wavelength currently used in the experiment is 532 nm. The two laser 
beams identical in divergence, intensity phase and polarization are obtained 
using a flat beam splitter. From experiment jlOj the value of the diffusion 
coefficient of the bead along the microtubule is obtained, D = 315nm^s~^. 

In this paper we study the heat signaling phenomena inside the micro- 
tubule excited by ultra-short laser pulses. 

2 Heat pulse transport on the molecular scale 

Molecular electronics is a new field of science and technology, which is evolv- 
ing from the convergence of ideas from chemistry, physics, biology, electronics 
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and information technology. It considers, on the one hand molecular mate- 
rials for electronic/optoelectronic applications, on the other hand attempts 
to build electronics with molecules at the molecular level. It is this second 
viewpoint which concerns us here: describing the heat transport at the level 
of single or few molecules. 

The heat and charge transport phenomena on the molecular scale are the 
quantum phenomena and electrons constitute the charge and heat carriers. 
We argue, that to describe the heat transport phenomena on the molecular 
level the quantum heat transport equation is a natural reference equation. 
The quantum heat transport equation for the atomic scale has the form 

where r is the relaxation time, D is the heat diffusion coefficient and T 
denotes temperature. For the atomic scale heat transport, the relaxation 
time equals 

^ 1 

Vh = —=ac, (2) 



where m is the electron mass and Vh is the velocity of the heat perturbation. 
Moreover, on the atomic level the temperature field T(r) is quantized by 
a quantum heat of energy, the heaton. The heaton energy equals 

Eh = rrievl. (3) 

Due to formula Q, the heaton energy is the interaction energy of electro- 
magnetic field with electrons (through the coupling constant a). 

At the molecular level we seek energy of interactions of the electromag- 
netic field with a molecule. This energy is described by the formula 

i^r = «^^mec^ (4) 

nip 

where me,mp are the masses of electron and proton, respectively. 

Considering the general formula for heaton energy (jH)), one obtains from 
formula (0]) for velocity of the thermal perturbation 

/ \ 1/2 
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Comparison of formulas (j2I) and (0) shows, that Vh scales with ratio (mg/mp)^/^ 
when the atomic scale is changed to the molecular scale; Vh is the Fermi ve- 
locity for molecular gas. 

Quantum heat transport equation ((H) has as a solution, for short time 
scale, (short in comparison to relaxation time r) the heat waves which prop- 
agate with velocity Vh- One can say that on the molecular level, the heat 
waves are slower in comparison to the atomic scale. 

From formulas (j21) and (0), the relaxation time can be calculated 

m„ h 

T = — 

It occurs, that relaxation time on the molecular scale is longer (ratio mp/me) 
than the atomic relaxation time. For standard values of the constants of the 
Nature 

a = me = 0.511 MeV/c^ rrip = 938 MeV/c^ (7) 

one obtains the following numerical values for Vh, t and Eh. Vh = 0.05 nm/fs, 
r = 44fs and = 10~^ eV. With those values of Vh and r, the mean free 
path 

A = TVh (8) 

can be calculated and A = 2.26 nm. It is interesting to observe, that in the 
structure of the biological cells (e.g. microtubules), some elements have the 
dimension of the order of the nanometer (e.g. microtubules). 

With the help of the heaton energy one can define the heaton tempera- 
ture, i.e., the characteristic temperature of the heat transport on the molec- 
ular scale, viz.: 

Tm = a^—rrieC^ 1.16 10^° K ^ 10~^a^meC^ ~ 316 K. (9) 
nip 

This defines what we generally term "room temperature" . At temperatures 
far below T^, the hydrogen bond becomes very rigid and the flexibility of 
atomic configurations is weakened. Most substances are liquid or solid below 
Tm- Biology occurs in environments with ambient temperature within an 
order or magnitude or so of Tm- 

In the following we assume that the ambient temperature is the initial 
temperature of the nerve cells. 

4 



(6) 



Figure 1: (a) The solution of QHT equation (0) for the following input pa- 
rameters Vh = 5 10-2 nm/fs, r = 44 fs, Tq = 300 K and At = 0.2r. (b) The so- 
lution of QPT |TT|) with the same input parameters. 



5 



Figure 2: (a, b) The same as in Fig. Ufa, b) but for At = r. 
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In Figs. HJ 121 the results of the theoretical calculations for the quantum 
heat transport on the molecular scale are presented. In Fig. ^a) the solution 
of Eq. (0) in one-dimensional case 

for the following input parameters Vh = 0.05 nm/fs, r = 44 fs, Tq = 300 K 
(initial temperature) and At duration of laser pulse = 0.2r is presented. 

In Fig. Wljo) the solution of quantum parabolic heat transport equation 
(QPT), (Fourier equation) 

dT d^T , , 

with the same input parameters is presented. 

In Figs. |21^a, b), the solutions of Eqs. ((T)), (fTT|l for the same input param- 
eter but for At = r are presented. 

From the analysis of the solutions of hyperbolic and parabolic quantum 
heat transport equation, the following conclusions can be drawn. In the 
case of QHT equation the thermal wave dominates the heat transport for 
At = 0.2 r. The finite value of Vh involves the delay time for the response of 
the molecular system on the initial temperature change. In the case of QPT, 
instantaneous diffusion heat transport is observed. From the technological 
point of view, the strong localization of thermal energy in the front of thermal 
wave is very important. 



3 Electron thermal relaxation in microtubules 

Clusters and aggregates of atoms in the nanometer size (currently called 
nanoparticles) are systems intermediate in several respects, between simple 
molecules and bulk materials and have been objects of intensive 
work. 

In this paper, we investigate the thermal relaxation phenomena in the 
nanoparticles - microtubules in the frame of quantum heat transport equa- 
tion. In the book [5 , the thermal inertia of materials heated with laser pulses 
faster than the characteristic relaxation time was investigated. It was shown, 
that in the case of the ultra-short laser pulses the hyperbolic heat conduction 
(HHC) must be used. For microtubules the diameters are of the order of the 



7 



electron de Broglie wave length. In that case to description of the transport 
phenomena quantum heat transport must be used [Hj. 



d^T dT /i 
ot-^ at m 

where T denotes the temperature of the heat carrier, r is the relaxation time 
and m denotes the mass of heat carrier. The relaxation time r is defined 
as [S] 

h 

r = 13 
where Vh is the thermal pulse propagation velocity 

Vh = (14) 

In formula ()14j) a is the coupling constant (for electromagnetic interaction 
a = /he) and c denotes the light velocity in vacuum. Both parameters r 
and Vh completely characterize the thermal energy transport on the atomic 
scale and can be named as atomic relaxation time" and "atomie" heat ve- 
locity. 

Both r and Vh are build up from constant of Nature, a, c. Moreover, on 
the atomic scale there is no shorter time period than r and smaller velocity 
build from constant of the Nature. In consequence, one can name r and Vh as 
elementary relaxation time and elementary velocity, which characterize heat 
transport in the elementary building block of matter, the atom. 

In the following, starting with elementary r and Vh, we intend to describe 
thermal relaxation processes in microtubules which consist of the parts 
(molecules) each with elementary r and Vh- To that aim, we use the Pauli- 
Heisenberg inequality 

ArAp > N^h, (15) 

where r denotes characteristic dimension of the nanoparticle and p is the 
momentum of energy carriers. The Pauli-Heisenberg inequality expresses the 
basic property of the iV-fermionic system. In fact, compared to the standard 
Heisenberg inequality 

ArAp > n, (16) 

we notice that, in this case the presence of the large number of identical 
fermions forces the system either to become spatially more extended for 
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fixed typical momentum dispension, or to increase its typical momentum 
dispension for a fixed typical spatial extension. We could also say that for 
a fermionic system in its ground state, the average energy per particle in- 
creases with the density of the system. 

A picturesque way of interpreting the Pauli-Heisenberg inequality is to 
compare Eq. (jl5|) with Eq. and to think of the quantity on the right 
hand side of it as the effective fermionic Planck constant 

hf{N) = Nh. (17) 

We could also say that antisymmetrization, which typifies fermionic am- 
plitudes amplifies those quantum effects which are affected by Heisenberg 
inequality. 

According to formula (fT7j) . we recalculate the relaxation time r, for- 
mula (fT^ and thermal velocity Vh, formula ()14|) for nanoparticle consisting 
fermions 

h-^hf{N) = Nh (18) 



and obtain 



r' = -7Tr,='^^- (20) 

Number A^ parts in nanoparticle (sphere with radius r) can be calculated 
according to the formula (we assume that density of nanoparticle does not 
differ too much from bulk material) 

N = (21) 

and for spherical with semiaxes a,b,c 

^abcpAZ , , 

N = ^ , (22) 

where p is the density of the nanoparticle, A is the Avogardo number, p is 
the molecular mass of particles in grams and Z is the number of the valence 
electrons. 
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With formulas (fTUI) and ((201), calculated de Broglie wave length X^^ 
and mean free path Xl^jp for nanoparticles 

\i = -^ = NhB, (23) 
>^Lfr = ^Lrt{ = Nh^f„ (24) 

where Xb and Xmfp denote the de Broglie wave length and mean free path 
for heat carriers in nanoparticles. 



4 Quantum transport in microtubules 

Microtubules are essential to cell functions. In neurons microtubules help 
and regulate synaptic activity responsible for learning and cognitive func- 
tions. While microtubules have traditionally been considered to be purely 
structural elements, recent evidence has revelated that mechanical, chemical 
and electrical signaling and communication function also exist, the result of 
microtubule interaction with membrane structures by linking proteins, ions 
and voltage fields respectively. 

The characteristic dimensions of the microtubules, crystalline cylinder 
10 nm in inner diameter are of the order of the de Broglie length for electrons 
in atoms. 

When the characteristic length of the structure is of the order of the de 
Broglie wave length, then the signaling phenomena must be described within 
the quantum transport theory. 

For quantum transport phenomena in microtubules we will apply the 
equation with the relaxation time described by formula (j2I) 

_ 2h _ h 
mf 2 E 

The relaxation time is the decoherence time, i.e. time until collapse of the 
wave function occurs, when the transition classical quantum phenomena 
is considered. 

In the following we will consider the time r for the atomic and multiatomic 
phenomena. As was shown in section 2 for atomic phenomena 

Ta ~ 10"^' S (25) 



10 



and when we consider multiatomic transport phenomena, with equal num- 
ber of agregates involved in transport we have, formula (j2Uj) 

TN = NTa. (26) 

The Penrose-Hameroff Orchestrated Objective Reduction (Orch OR) model [2] 
proposes that quantum superposition - computation occurs in nanotubule au- 
tomata within brain neurons and glia. Tubulin subumits within microtubules 
act as qubits, switching between states on a nanosecond scale, governed by 
London forces in hydrophobic pockets. These oscillations are tuned and or- 
chestrated by microtubule associated proteins (MAPs) providing a feedback 
loop between the biological system and the quantum state. These qubits 
interact computationally by nonlocal quantum entanglement, according to 
the Schrodinger equation with preconscious processing continuing until the 
threshold for objective reduction (OR) is reacted {E = ^). At that instant, 
collapse occurs, triggering a "moment of awareness" or a conscious event - 
an event that determines particular configurations of Planck scale experi- 
ential geometry and corresponding classical states of nanotubules automata 
that regulate synaptic and other neural functions. A sequence of such events 
could provide a forward flow of subjective time and stream of consciousness. 
Quantum states in nanotubules may link to those in nanotubules in other 
neurons and glia by tunnelling through gap functions, permitting extension 
of the quantum state through significant volumes of the brain. 



Table 1: 



Event 


T [ms] 


E 


N 

number of aggregates 


T [ms] 


Buddhist moment 


13 


4 ■ 10^5 nucleons 


10^5 


10 


of awareness 










Coherent 40 Hz 


25 


2 ■ 10^5 




10 


oscillations 










EEG alpha rhytm 


100 




10^4 


1 


(8 to 12 Hz) 










Libet's sensory 


500 






1 


threshold 











From E = j;, the size and extension of Orch OR events that correlate with 
subjective or neurophysiological description of conscious events can be calcu- 
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lated. In Table 1 the calculated T (Penrose-Hameroff) and r— formula (|26|) 
are presented 0- 



5 Heisenberg uncertainty principle for ther- 
mal phenomena in microtubules 

Efficient conversion of electromagnetic energy to particle energy is of funda- 
mental importance in many areas of physics. The nature of intense, short 
pulse laser interactions with single atoms and solid targets has been subject 
of extensive experimental and theoretical investigation over the last 15 years. 
Recently, the interaction of femtosecond laser pulses with Xe clusters was 
investigated and strong X-ray emission and multi-keV electron generation 
were observed. Such experiments have become possible, owing to recently 
developed high peak power lasers which are based on chirped pulse am- 
plification and are capable of producing focused light intensity of up to 
1014 _ 1019 Wcm-2. 

In intensely irradiated clusters e.g. microtubules, optically and collision- 
ally ionized electrons undergo rapid collisional heating for short time (< Ips) 
before the cluster disintegrates in the laser field. Charge separation of the 
hot electrons inevitably leads to a very fast expansion of the cluster ions. 
Both electrons and ions ultimately reach a velocity given by the speed of 
sound of the cluster plasma. 

When the intense laser pulse interacts with atomic clusters ionization to 
very high charge states is observed. The high Coulomb field certainly influ- 
ences the thermal processes in clusters. In the chapter, the new QHT equa- 
tion is formulated in which the external — not only thermal forces are in- 
cluded. The solution of the new QHT for Cauchy boundary conditions will 
be derived. The condition for the distortionless propagations of the thermal 
wave will be formulated. 

Now, we develop the generalized quantum heat transport equation which 
includes the potential term. In this way, we use the analogy between the 
Schrodinger equation and quantum heat transport equations. Let us con- 
sider, for the moment, the parabolic heat transport equation with the second 
derivative term omitted 

— = -V^T. (27) 

ot m 

When the real time t —>■ ^ and T ^ \E', Eq. (j27p has the form of a free 
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Schrodinger equation 



dt 2m ^ ^ 



The complete Schrodinger equation has the form 

where V denotes the potential energy. When we go back to real time t — > —2it 
and —y T, the new parabolic quantum heat transport is obtained 

or 2V^ , , 

= — V^T - — T. 30 
at m a 

Equation (jHUj) describes the quantum heat transport for At > r. For heat 
transport initiated by ultra-short laser pulses, when At < r one obtains the 
generalized quantum hyperbolic heat transport equation 

d^T dT h^,^ 2V ^ , , 

^arH-g7 = -V=r--T. (31) 

Considering that r = h/mv"^, Eq. (jHT|) can be written as follows: 

1 d'^T mdT 2Vm^ , , 

Equation ()32j] describes the heat flow when apart from the temperature gra- 
dient, the potential energy V operates. 

In the following, we consider the one- dimensional heat transfer phenom- 
ena, i.e. 

I d'^T mdT 2Vm„ d'^T 



dt'^ ~^ h dt ~^ h dx^ ^"^"^^ 

For quantum heat transfer equation we seek solution in the form 

T(x,t) = e-*/=^"M(x,t). (34) 

After substitution of Eq. (jMj) into Eq. (jH^ . one obtains 

1 d^u d'^u , , ^ 

where 

2'Vm /mv 
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In the following, we will consider the constant potential energy V = Vq. 
The general solution of Eq. for Cauchy boundary conditions, 



.(x,0) = /(x), 



dt 

has the form [Hj 



= Fix), (37) 

i=0 



, , f(x-vt) + f(x+vt) 1 r^''* , 

u{x, t) = ^ ^ + - / $(a;, t, z)dz, (38) 



where 



<f(x,t,z) = -F{z)Jo (-y/{z - X)2 - + 6t/(^)_^ 



{z — x)^ — vH"^ 



2h / ^2 



t;^ (39) 



and Jq{z) denotes the Bessel function of the first kind. Considering formu- 
las (j3l|l . (j35|) . (j36|l the solution of Eq. (j33|l describes the propagation of the 
distorted thermal quantum waves with characteristic lines x = ±vt. We can 
define the distortionless thermal wave as the wave which preserves the shape 
in the field of the potential energy Vq. The condition for conserving the shape 
can be formulated as 

2Vm /mf\2 



-(ZY- (40) 



\2hJ 

When Eq. holds, Eq. ^ has the form 

dt^ dx^' ^ ^ 

Equation (PT|) is the quantum wave equation with the solution (for Cauchy 
boundary conditions (jSZI) 

, , f(x-vt) + f(x + vt) 1 , , , 

u{x, t) = '-^ '—11 L + / F{z)dz. (42) 

^ Jx-Vt 

It is quite interesting to observe, that condition pn]) has an analog in the 
classical theory of the electrical transmission line. In the context of the 
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transmission of an electromagnetic field, the condition q = describes the 
Heaviside distortionless line. Eq. (^01) — the distortionless condition — can 
be written as 

Vqt ~ h, (43) 

We can conclude, that in the presence of the potential energy Vq one can 
observe the undisturbed quantum thermal wave only when the Heisenberg 
uncertainty relation for thermal processes ()43|) is fulfilled. 

The generalized quantum heat transport equation (GQHT) (jH^ leads 
to the generalized Schrodinger equation. After the substitution t — ^ it/2, 
T — i> in Eq. (jHSl), one obtains the generahzed Schrodinger equation (GSE) 

Considering that r = h/mv'^ = h/ma^c^ {a = 1/137 is the fine-structure 
constant for electromagnetic interactions) Eq. can be written as 

in— - -—V^xiJ + - (45) 

dt 2m ma^c^ df^ 

One can conclude, that for time period At < h/ma^c^ ~ 10^^'' s the de- 
scription of quantum phenomena needs some revision. On the other hand, 
for At > 10^^^ in GSE the second derivative term can be omitted and as the 
result the SE is obtained, i.e. 

ih— = V^^ + V^. (46) 

dt 2m ^ ^ 

It is quite interesting to observe, that GSE was discussed also in the context 
of the sub-quantal phenomena. 

Concluding, a study of the interactions of the attosecond laser pulses with 
matter can shed the light on the applicability of the SE to the study of the 
ultra-short sub-quantal phenomena. 

The structure of the Eq. ()35|) depends on the sign of the parameter q. 
For quantum heat transport phenomena with electrons as the heat carriers 
parameter q is the function of potential barrier height (Vq) and velocity v. 
Considering that velocity v equals 

1 nm , , 

w = ^ac=1.26— , 47 
V3 fs 
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Figure 3: Parameter q (formula ()36|)) as the function of the barrier 
height (eV). 



parameter q can be calculated for typical barrier height Vq > 0. In Fig. El 
the parameter q as the function of Vq is calculated. For g < 0, i.e., when 
Vq <1.25 eV, Eq. (jH^ is the modified Klein-Gordon equation. 
For Cauchy initial condition 



u{x,o) = /(x), 



du{x, a) 



dt 



the solution of the Eq. (j35p has the form 

f{x-vt) + f{x + vt) 



u{x, t) 



1 

+ 2^ 



+ 



2 

x+vt 
x~vt 



q)t 



x+vt 

/(C)- 



^/-qivH^-ix-cy) 



(48) 



(49) 



x—vt 



When g > Eq. ()35p is the Klein-Gordon equation (K-G) well known from 
application to elementary particle and nuclear physics. 



16 



For Cauchy initial condition (j48j) . the solution of (K-G) equation can be 
written as 



u{x, t) 



f{x-vt) + f{x + vt) 



2 

-| px+vt 



(50) 



x+vt 



/(C)- 



4 



^q{vH^ - (x - CY) 



x—vt 



^vH'^ -{x- C)^ 



Both solutions (jl^ and (jMIj) exhibit the domains of dependence and influ- 
ence on modified Klein- Gordon equation and Klein-Gordon equation. These 
domains, which characterize the maximum speed at which thermal distur- 
bance travel are determined by the principal parts of the given equation (i.e., 
the second derivative terms) and do not depend on the lower order terms. 
It can be concluded that these equation and the wave equation (for m = 0) 
have identical domains of dependence and influence. 



6 Thermal wave packets induced by attosec- 
ond laser pulses 

Equation (|35|) is the modified Klein-Gordon equation which can be written 
as 

Du-(^yu{x,t) = 0, (51) 

where d'Alembert operator □ is equal 

1_ d d 



□ 



V dt^ dx^ 



(52) 



The ordinary Klein-Gordon equation for the particle with mass m and ve- 
locity V is of the form 



(53) 



Equation ()51|) can be split into its real and imaginary parts. Putting for 
u{x^ t) 



u{x, t) = 9ft(t, x) exp 



h 



S{t,x) 
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one obtains 

V^\d^S){d,S) = - , (54) 

where 

r]ab = diag(l, -1, -1, -1), a,b=l, 2, 3, 4. 

We use Mackinnon's suggestion [H] therefore we look for solutions that 
satisfy the equation 

□3? r mv " 



In this way Eq. ()54|) becomes 



(55) 



V'''{daS){dbS) = = P^P>^. (56) 

If the velocity v is constant, we have from (jMj) 

S = -P^P^ (57) 

and the de Broglie relation P'" = HK'" hold where is the wave num- 
ber and P^ is the classical relativistic four momentum, P^ = {^p}- In 
paper [6] L. Mackinnon constructs a wave packet considering that a wave 
$ = $0 exp [iut] of frequency u = is associated with a particle of rest 
mass m and that for an observer moving with a constant velocity v with 
respect to the particle, the associated wave (be means of a Lorentz transfor- 
mation) acquires the form $ = $oexp ^iu'j (t — ^a;)]. Mackinnon showed 
the wave packet is compatible with the basic experiments of quantum me- 
chanics and does not spread in time the found 

5J = A^, (58) 

gr 

where A = constant and g = and g = and 

r = 7(0; — vt) (59) 
is the distance from the particle portion, so that 

(60) 



, , . sm gr 

u{x,tj = A exp 

gr 



- — (i?t — px) 
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is the Lorentz boost of the solution 



$ = A — exp \-tujt ] , 

gr' 

where r' = x. This solution was first found by de Broglie [S] and also used in 
the stochastic interpretation of quantum mechanics by Vigier and Gueret P . 



19 



References 

[1] J. G. Lambert et al., Cell. Mol. Biol, 45, (1999), p. 905. 

[2] R. D. Vale et al., Cell, (1985), p. 449. 

[3] D. G. Weiss et al.. Cell Motil. Cytoskeleton, 6, (1986), p. 128. 

[4] D. A. Smith, R. M. Simmons, Biophysical Journal, 80, (2001), p. 46. 

[5] M. Kozlowski, J. Marciak-Kozlowska, From Quarks to Bulk Matter, 
Hadronic Press, 2001. 

[6] S. R. Hameroff, J. Tuszynski, Conscious events as orchestrated space- 
time selections. Journal of Consciousness Studies, 3, (1996), p. 36. 

[7] L. Mackinnon, Letter at Nuovo Cimento, 32, (1981), p. 311. 

[8] L. de Broglie, C. R. Acad. Sci., 180, (1925), p. 498. 

[9] J. R Vigier and Ph. Gueret, Found. Phys., 12, (1982), p. 1057. 

[10] G. Cappclo et al., Phys. Rev., E68, (2003), 021907-1. 



20 



MTE[q<0], WE[q=0], K-GE[q>0] 




Barrier height [eV] 



(a) QH-TEMPERATURE T0=300 K 
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